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Abstract 

We show that hysteresis can occur in the transition between a neutral plasma 
and the insulating gas consisting of neutral pairs bound by Coulomb attrac- 
tion. Since the transition depends sensitively on the screening length in the 
plasma, regions of bistability occur in density-temperature phase space. We 
present numerical results which indicate where these regions occur for sys- 
tems such as spin-polarized hydrogen, positronium gas, and excitons in a 
semiconductor. 
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The system consisting of eqnal nnmbers of positively and negatively charged spin-1/2 
particles provides a snrprisingly rich snbject for the stndy of phase transitions. At high 
density, the system is the well-known nentral Fermi plasma. As density decreases, the 
system becomes nnstable to the formation of bonnd complexes which behave as bosons. Un- 
der certain conditions at low temperatnre, these qnasi-bosons can form a Bose condensate. 
Presently, experimental searches are nnder way or proposed for observation of Bose conden- 
sation in several systems of this type, namely excitons in semicondnctors |l|], hydrogen 
and positroninm [Q . Depending on the mass ratio and spin degeneracies of the positive and 
negative particles, other phases snch as solid and Fermi liqnid can arise. 

The phase transition from condncting plasma to insnlating bonnd-pair gas, known as the 
Mott transition, has been discnssed in the literatnre for well over three decades. The general 
condition for this transition is that the screening length in the plasma be comparable to or 
less than the intrinsic bonnd-state Bohr radius. A point that has been less well appreciated 
is the possibility for bistability and hysteresis in this transition. The basic reason for this 
hysteresis is that neutral bound states will not contribute to long-range screening. Therefore, 
if the gas is initially in the insulating phase, a transition to conducting plasma will occur 
only when the number of free charged particles created by thermal dissociation of bound pairs 
becomes high enough that the screening due to these particles causes bound states to become 
unstable. In general, the density of this transition, called the “ionization catastrophe,” 
is not the same as the Mott transition in the reverse direction that occurs when the system 
is initially in the plasma state. 

In this paper, we present a simple model that illustrates the hysteresis in this kind of 
system. Much of the discussion is in terms of equal-mass positive and negative particles 
(excitons or positroninm) but the discussion can be generalized to apply to the large-mass- 
ratio (hydrogen atom) case as well. For the sake of discussion, the bound states will be called 
excitons. In all of the following we assume nondegenerate (Maxwell-Boltzmann) statistics 
for both bound states and free particles. This will limit the regions of phase space where 
the model is valid (i.e. it will not be valid near the Bose-Einstein phase transition or for the 
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degenerate Fermi gas), but the primary effect of hysteresis does not depend on the statistics. 
We also ignore van der Waals interactions between bound states and the formation of bound- 
state complexes. 

Rate equations. Instead of the usual partition function Z = “Planck- 

Larkin” partition function |§,^, Z = _j_ A/fc^T), is used, which takes into 

account the fact that states with A -C ksT are quasi-free. This first-order correction to 
the Boltzmannian partition function was initially calculated in order to keep the density of 
states of atoms near the n ^ oo continuum from an unphysical infinity. In the context of 
the neutral plasma system studied here, it enforces the physical result that when the binding 
energy is zero, the population in bound states must become identically zero. Ebeling et ah 
P have used first-order perturbation theory to show that this partition function has sound 
physical basis. 

For species conversion via three-body collisions, the rate of change of the bound state 
(exciton) population n^x is then given as follows: 
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where iie is the free electron density, with Uex + rie = n, the total pair density; Qe, Qh and 
gex are the electron, ion (hole) and exciton spin degeneracy, respectively, and m and fi are 
the total and reduced mass, respectively. The “quantum density of states” is given by 
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where A/? is the deBroglie wavelength. Setting dnex/dt = 0 in (P with A ;§> ksT gives 
the standard Saha equation {nl/nq) = {geQh/ gex){lJ^'/ . The constant A, 

which depends on the cross section and may be temperature dependent, determines the 
absolute rate of conversion and does not enter into steady-state calculations. 

The interesting dynamics in this system come from the dependence of A on ng. The 
binding energy (Rydberg) of a single bound state (exciton) is given by Aq = h^/2a^/i, 
with the Bohr radius a = h^e/e^g,. (The dielectric constant e is included in the case of a 
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polarizable medium, e.g. for excitons in a semiconductor.) When screening is present, so 
that the particles do not interact with pure Coulomb attraction but with a Yukawa-type 
interaction, the binding energy decreases approximately according to p| 



A(ne,T) = 



Ao (l - t+(|ipt) < 1 
0 qa >1. 



(3) 



For 1/q ~ a, variational calculations give slight (20%) corrections to the above formula 0 
which are not of interest here. The screening constant is given by the Debye-Huckel formula, 

m 
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It is important to recognize that only free charged particles contribute to screening- the 
bound states are neutral. 

Numerical results. It is convenient to set Ue = ny'^ and use fiex = — hg to write the 
dynamics in dimensionless form 
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dr f(y) L 

where r = AnH/2 is a scaled time, and f{y) = exp{F{y)) — {1 + F{y)) is the partition 
function in terms of 

{9e9h/gf f 
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The parameters a and f3 determine the temperature and density 
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There is an equilibrium at j/ = 0 which is unstable due to thermal ionization of excitons. 
In addition, there is always at least one stable equilibrium described by 



0 = 



«(1 - 2 /^) - 2 /^/( 2 /) 



(9) 
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For appropriate choices of n and T there can be three solutions; two stable and one unstable. 
The locus of points (n, T) dehning the region of bistability has been determined numerically. 
When there is a fully ionized plasma state at y = 1 then F{1) = 0, and the model fails to 
quantitatively describe the approach to y = 1 due to the divergence of as F — > 0. 

Figure 1 shows numerical results when = 2 and gex = 1- The unitless pair 

density na^ and the unitless temperature {y/m){a^/Xj^) are control parameters; only the 
fraction of free carriers y'^ = nejn can vary. The binding energy A(ne) is calculated using 
a self-consistent value of Ue- The solid line marks the boundary between the region with 
one steady-state solution and the region with three steady-state solutions, one of which is 
unstable. 

Figure 2 illustrates the behavior of the system in these regions, by plotting the ionized 
fraction as a function of total pair density at hxed temperature. For temperature well above 
the three-solution region (Fig. 2(a)), the system remains mostly ionized at all densities. At 
lower temperatures (Figs. 2(b) and (c)), a sharp transition appears between a pure plasma 
and exciton-rich phase. Finally, when temperature is lowered further, a classic hysteresis 
curve occurs (Fig. 2(d)). Two different solutions for the ionized fraction are stable against 
fluctuations. 

Regions of phase space. By simple arguments, we can see how hysteresis occurs in this 
system. Figure 3 shows the relevant crossovers for the case m = 4/i and ge = 9h = 2, gex = 1 
of Figs. 1 and 2. If we assume that the system is in the plasma state, then we expect a Mott 
transition when 1/a = q where q is given by the Debye formula (||) with Ug = n. In terms of 
the unitless parameters na^ (unitless pair density) and a^/A|, (unitless temperature), this 
becomes 

^ ( 10 ) 

shown as the heavy solid line in Figure 3. On the other hand, if we assume the system is in 
the insulating gas (excitonic) phase, we would expect a transition when 1/a = q where q is 
given by (^), but instead of setting Ug equal to n, we would use the equilibrium value from 



5 



(I), assuming A = Aq and Ue -C n. This yields 
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shown as the light solid line in Figure 3. Since at very low temperatures the number of 
thermally ionized particles becomes extremely small, the “ionization catastrophe” density 
becomes exponentially large at low temperature. For comparison, the conditions for a tran- 
sition to a Fermi degenerate gas {na^ ~ 1) and for a transition to a weakly-interacting Bose 
condensate {na^ = 2.612(a^/A|))^/^) are shown as the dashed and dotted lines in Figure 1, 
respectively. Since the model presented here involves only Maxwellian statistics, it is only 
valid for densities well below both of these curves. The presence of the plasma state will 
prevent Bose condensation, however, so the proximity of the Mott transition to the Bose con- 
densation boundary in certain regions of phase space merits attention in all composite-boson 
condensation experiments. 

At low temperature and low density, the system is an insulating excitonic (bound-state) 
gas, assuming repulsive interactions between excitons. As density is raised, in the low- 
temperature limit the system remains an insulating gas until an ionization catastrophe 
occurs. As seen in Fig. 3, this can occur at substantially higher density than the Mott 
transition given by Eq. (|l0|) . If the system starts in the plasma state in this temperature 
regime, however, it will remain a plasma indehnitely unless the density falls below the Mott 
density given by (pT|). 

These two curves essentially form the boundaries of the region of hysteresis shown in 
Fig. 1. (The curve in Fig. 3 given by (^) lies above the phase boundary for an ionization 
catastrophe found numerically in Fig. 1 because it does not use a self-consistent value 
of the binding energy A(ne).) At temperatures above the intersection of the solid curves, 
however, the system undergoes a single phase transition to a plasma state, since an ionization 
catastrophe occurs at lower density than the Mott transition given by (^). At very high 
temperature, the “ionization catastrophe” is not a clearly dehned transition, since even at 
low density a substantial fraction of the system is ionized. 
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Application to experiment. Preliminary results from experiments with excitons in 



the semiconductor CU 2 O seem consistent with this picture. As reported in Ref. |^, excitons 
with binding energy A = 150 meV exist in this material and have been observed at room 
temperature. In those experiments, the value of na^ was roughly 0.003 (pair density ~ 10^® 
cm“^, excitonic Bohr radius 7 A), a ? / was approximately 0.04 (T = 300K, mass m = Suiq), 
and /i/m ~ 4 (roughly equal electron and hole mass.) This puts the system right in the region 
of the crossover seen in Figures 1 and 3. In the experiments, excitons appeared only when 
created resonantly by laser excitation. When free carriers were created by laser excitation 
with photon energy well above the band gap, no exciton luminescence ever appeared. In 
other words, two different stable states occurred at roughly the same temperature and pair 
density. 

Clearly, there is room for more experiments and theory on this interesting system. Since 
the model presented here assumes spatial homogeneity, pattern formation and related non- 
linear dynamical effects do not occur. The possibility exists that when the effects of spatial 
variations are taken into account, regions of coexistence and very interesting chaotic behavior 
may be seen in the exciton/plasma phase transitions. 

These results also relate to experiments on phase transitions of atoms with attractive 



interactions, e.g. the search for a metallic phase of hydrogen at high pressure |]l3- Because 
of attractive interactions, the lower, insulating state may not be a gas, as assumed here, 
but instead a solid or liquid state. In many cases, however, screening is assumed to play 
an important role in a transition to a conducting phase. The results here indicate that the 
densities necessary for a spontaneous transition to such a phase may be much higher than 
expected from equilibrium calculations, since a metastable solid or liquid state comprised of 
bound pairs will have poor screening even at high density. 
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FIGURES 



FIG. 1. The boundary between the regions of phase space with one steady-state solution 
and three steady-state solntions of Eq. (D in the text, for the case ge = gh = ‘2, g = 1. 

FIG. 2. Steady-state solutions of n^ln from Eq. (^, as a fnnction of na^, for four 
temperatures. 

FIG. 3. The regions of phase space indicated by the simple equations ( 0 ) and (pdD in 
the text. The model breaks down above the critical density for Bose condensation (dotted 
lines, for two different values of m//x corresponding to excitons and hydrogen) and in the 
Fermi degenerate region {no? ~ 1, indicated by the dashed line.) 
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